Anomalously Soft Droplets and Aging in Short-ranged Spin Glasses 
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We extend the standard droplet scaling theory for isothermal aging in spin glasses assuming that 
the effective stiffness constant of droplets as large as extended defects is vanishingly small. A novel 
dynamical order parameter and the associated dynamical susceptibility xd emerge. Breaking of 
the fluctuation dissipation theorem takes place at qo well below the equilibrium Edwards- Anderson 
order parameter qea(> ?d) at which the time translational invariance is strongly broken. The 
scenario is examined numerically in a 4 dimensional EA Ising spin glass model. 



Recently aging in spin glasses ||, [| [|, || has at- 
tracted renewed interest both in experimental, numeri- 
cal and theoretical studies. The phenomena provide very 
interesting and almost the only realistic way to explore 
large scale, non-trivial properties of glassy phases. A 
major theoretical progress was the development of dy- 
namical mean-field theories (MFT) for spin glasses and 
related systems || which lead to findings of some novel 
view points in glassy dynamics such as the concept of 
effective temperature [Q. However, the MFT themselves 
do not provide insights into what will become important 
in realistic finite dimensional systems, such as nucleation 
processes. 

The droplet theory [||, ^, [l(| developed a scaling the- 
ory assuming droplet excitations as basic nucleation pro- 
cesses. It provides a perspective for equilibrium and 
dynamical properties of finite dimensional spin glasses. 
Though it was developed more than a decade ago, many 
of its predictions have remained to be clarified. 

Recently, anomalous low energy and large scale excita- 
tions were found [ pd| in active studies of spin-glass mod- 
els at T = 0. In the present letter, we develop a refined 
scenario for isothermal aging in spin glasses where such 
anomalous low energy excitations play a significant role. 
Most importantly, it explains the fundamental experi- 
mental observation that the field cooled (FC) suscepti- 
bility is larger than the zero field cooled (ZFC) suscep- 
tibility [|[ [|. Such an anomaly was reproduced success- 
fully in the dynamical MFT but not in the conventional 
droplet theory. To simplify notations, we discuss sys- 
tems of N Ising spins Si = ±1 (i = 1, . . . , N) in a d- 
dimensional space coupled by short-ranged interactions 
of energy scale J with random signs with no ferromag- 
netic or anti-ferromagnetic bias. Extension to other types 
of spin glasses can be done straightforwardly. 

Fisher and Huse |l(| noticed that droplet excitations 
can be softened in the presence of a frozen-in defect or 
domain wall compared with ideal equilibrium with no 
extended defects. This is because droplets which touch 
the defects can reduce the excitation gap compared with 
those in equilibrium. Let us consider system with a 



frozen-in defect of size R: a large droplet of size R is 
flipped with respect to T, which is a ground state of an 
infinite system, and then it is frozen. The typical free- 
energy gap F^ji of a smaller droplet of size L in the 
interior of the frozen-in defect is expected to scale as 



ptyp 
L,R 



T eS [L/R](L/L ) 



L < R 



(1) 



with Lq being a certain microscopic length scale and 
Y c ft[L/R] being an effective stiffness which is only 
a function of the ratio y = L/R. We assume as 
usual that the probability distribution of the free-energy 
gap Fl,r follows the scaling form p(Fi, t R)dFi„R — 
P{Fl,r/ 'F^^dFL^Rj F^x with non- vanishing amplitude 
at the origin p(0) > which allows marginal droplets. 
For y <C 1, T e g[y] will decrease with y as |1C[ 



T e fi[y]/T = 1 - c v y 



6,-6 



for 



y < 1. 



(2) 



Here T is the original stiffness constant T = T c ff(0). A 
basic conjecture, on which our new scenario based, is that 
at the other limit y ~ 1 the effective stiffness vanishes as, 



T eff [y]/T~(l-y)< 



(3) 



with < a < 1 being an unknown exponent, and that 
the lower bound for -F^ls should be of order J, say Fq. 

In spin glasses of finite sizes R, it is likely that the 
existence of boundaries will intrinsically induce certain 
defects as compared with infinite systems |l^, |l3|, pd[ . 
Then droplet excitations as large as the system size itself 
L — R may be anomalously soft as conjectured above. 
Such an anomaly is indeed found in recent studies [ jjT| . 
It explains the apparently non-trivial overlap distribu- 
tion function P(q) found in numerous numerical studies 
of finite size systems [ fl5|. A lthough there new exponent 
6' = was conjectured [|ll|, we consider it better to at- 
tribute this to the zero stiffness constant T e ff[l] = as 
in (Q). The stiffness exponent > 0, on the other hand, 
is associated with a defect in T as we adopt in (Q). 

To explain the consequence of our conjecture above in- 
troduced, let us consider thermal fluctuations and mag- 
netic linear responses of droplet excitations of size L 
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within the frozen-in defect at scale R. Each droplet 
excitation will induce a random change of the magne- 
tization of order Ml ~ my (L/L ) d where m is the 
average magnetic moment within a volume of L . The 
thermal fluctuation can be measured by an order param- 
eter q — N^^-^iiSi) 2 where the sum runs over sites 
in the interior of the frozen-in defect which contains 
Nd oc (R/L ) d spins. The magnetic linear response by 
weak external magnetic field h is measured by a linear 
susceptibility x = N^ 1 Y,t{ s i)h/h where J2i( S i)h is the 
induced magnetization by the field. In the absence of 
any droplet excitations, q = 1 holds due to the normal- 
ization of spins. The reduction from 1 due to a droplet 
excitation at scale L is of order Mf j p(0)(k B T / F^^) with 
k B being the Boltzmann constant. Correspondingly, the 
induced magnetization of droplets by h at scale L is of 
order M L p(0)(hM L / 'F^). 

Let us construct a toy droplet model || defined on 
logarithmically separated shells of length scales L/Lq = 
a k < R/Lq with a > 1 and < k < n 2 < ri\ where 
R/Lq = a" 1 and L m /Lo — a™ 2 . For each shell an optimal 
droplet is assigned whose free-energy gap is minimized 
within the shell. Then F^j, will be of order Fq at the 

shell k — n± so that we have F^^ = Fl^rO- — Sk >ni ) + 
Fo5k,m ■ For simplicity, droplets at different scales are 
assumed to be independent from each other. Replacing 
the sum YlkLo by an integral J L m dL/L we obtain, 



l-q(L m ,R)=k B Tx(L m ,R) 
I- A a (ln(L/i?)) 



dL 
T 
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L,R 



= p(0)m 2 k B T x 

A„(ln(£/fl)) 
F 



(4) 



where A a (z) is a pseudo S- function of width In a [ ^7[ . 
Note that the fluctuation dissipation theorem (FDT) is 
satisfied. 

It is useful to consider assymptotic behaviour at large 
sizes R/Lq 3> I with the ratio x = L m /R being fixed. 
We obtain for < x < 1, 



q(xR,R) = q E A- 



p{0)m 2 k B T 
' T(R/Lo) e 



p{0)m 2 k B T 



A{x)- ^"l e a (ln») 
(5) 

with Q a (z) being a pseudo step- function of width In a 
and A(x) = Cdyy-^ 9 -f*dyy-i- e (r/r cS {y](l- 
A a (lny)) — 1). Note that the second integral converges 
because of (||) and the inequality < (rf-l)/2 Here 
qea is the usual Edwards- Anderson (EA) order parame- 
ter evaluated in the absence of any extended defects as 



lim lim q(L m , R) 

L m — >oo R — >oc 



1 — cp(0)m 



k B T 



2 K B 



(6) 



with c = Jj 00 dyy~ x ~ e . The order of the limits is crucial. 
Thus as far as < x < I, the order parameter converges 
to the equilibrium EA order parameter <7ea in the large 
size limit. The associated equilibrium susceptibility xea 
is defined as fceT^EA = 1 — Qea- 



In the intriguing case x ~ 1, A(x) will remain finite 
as far as < a < 1. At x ~ 1 the last term of (|^) 
contributes and we obtain, 



q(R, R) = q B + p(0)m 2 



k B T 



T(R/L 



(T) 



where we have defined the dynamical order parameter 

2 k B T 



q D = lim q(R, R) = <?ea - p(0)m 



(8) 



This is one of the main results of the present work. Nat- 
urally, we can define the associated dynamical linear sus- 
ceptibility xd as k B T\ = 1 — Id- As we discuss below qc 
and xd play significantly important role in the dynamical 
observables of aging. 
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FIG. 1: Different asymptotic regimes of the two-time quanti- 
ties. Here the asymptotic limit L(t w ) — > oo is considered with 
the ratio x = L(T)/L(t w ) being fixed. 

We are ready to consider a simple isothermal aging pro- 
tocol: the temperature is quenched [|l6| at time t = to a 
temperature T below its transition temperature T c from 
above T c and the relaxational dynamics i.e., aging is mon- 
itored. The dynamical magnetic susceptibility is gener- 
ally defined as x(Mw) = N- l (M(t))/h where (M(t)) is 
the total magnetization measured at time t and field h. 
In the thermoremanent magnetization (TRM) measure- 
ment the field is switched on for a waiting time i w and 
then cut off: h(t) = hd(t w — t). On the other hand, in the 
ZFC magnetization measurement, the field is switched 
on after the waiting time: h(t) = h0{t — i w ). We denote 
the susceptibilities of these protocols as xtrm (t, i w ) and 
Xzfc(^^w)- Another important quantity is the magnetic 
autocorrelation function, C{t,t w ) = N' 1 {M(t)M(t w )) = 
N~ 1 J2i(^i('t)^i('tv))- The last equation follows in the 
absence of ferromagnetic or anti-ferromagnetic bias. 

Within the droplet picture, aging proceeds by coarsen- 
ing of domain walls jl0| as usual phase ordering processes 
p8[ . Following [|To), the whole process may be divided 
into epochs such that the typical size of the domain is 
Lq, aLo, a 2 Lo, . . . , a n Lo, .... At each epoch, droplets of 
various sizes up to that of the domain can be thermally 
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activated or polarized by the magnetic field. Importantly 
the domain wall serves as the frozen-in extended defect 
for the droplets and reduces their stiffness constant. 

In the droplet picture [jlo| all dynamical processes, 
including both the domain growth and droplet excita- 
tions in the interior of the domains, are thermally ac- 
tivated processes which yield a logarithmic growth law 
of the time dependent length scale L(i). In practice, 
empirical power laws with temperature dependent ex- 
ponent work as wellju], However, the logarithmic 
growth law is supported by a recent experiment ]l9[] 
where the effects of critical fluctuations are considered 
(see also||c], ^l], ^3|). In the following, we focus on 
scaling properties of the two time quantities as functions 
of the dynamical length scale rather than on the growth 
law itself. 

As far as L(t) < L(t w ), the autocorrelation func- 
tion and the ZFC susceptibility are related to the gen- 
eralized order parameter and linear susceptibility de- 
fined in (|J) as C(r + t w ,t w ) — q(L(r), L(t w )) and 
kBTxzFc{T + t w ,t w ) = 1 - C(r + t w ,t w ) satisfying 
the FDT. Especially, the equilibrium correlation and re- 
sponse are obtained by taking the limit L(t w ) —> oo 
with fixed L(f). In this special limit one finds Jl0| , 
Ccq(r) = q EA + cp(O)m 2 (k B T/T)(L(T)/L o )' and 



In the aging regime L(t) ~ Lit) > L(t w ) we expect, 



Xeq(r) = lim xzfc(t + £ 



,*w) = Xea-c- 



p(0)m 2 



r(L(r)/L y 

(9) 

In the quasi-equilibrium regime L(r) < L(t w ) waiting 
time dependence or violation of TTI (Time translational 
invariance) is present in a weak manner as correction 
terms to the ideal equilibrium behavior (^J). The scaling 
form of the correction terms can be found using ( p[) w hich 
is relevant for relaxation of AC susceptibilities JlO| , [l9[ . 
Actually it allows one to set up numerical extrapolations 
to obtain the ideal equilibrium behavior [H], [2^]. In 
the limit L(t w ) — > oo with x = L(r)/L(t w ) < 1 being 
flexed, we find C — > q E A and xzfc — » Xea (see Fig. [l]). 
Much stronger violation of TTI shows up below qea- 

In the crossover regime L(t) ~ L(t w ), the anoma- 
lously soft droplets as large as the size of the domain 
L(t w ) are also needed to be taken into account. From 
(0) we immediately find, C(t +t w , iw)|z( T )~z(t w ) ~ <7d + 
A(l)p(0)m 2 (k B T/T)(L(t w )/L o )~ e where q B is the novel 
dynamical order parameter we introduced above. Thus 
as a function of x, there should be a vertical drop from 
qea to gD at x ~ 1 in the asymptotic limit L(t w ) — > oo. 
Correspondingly, the ZFC susceptibility should jump up 
vertically from xea to xd at x ~ 1 in the same asymp- 
totic limit (see Fig. |l|). Such an abruptness will be 
absent as functions of r/t w . Indeed it is well known 
in experiments that the so called relaxation rate 
S(t) = cIxzfc(t + t w ,t w )/dlog(T/T ) has a pronounced 
peak at around r ~ t w . 



C(t,t w )~ qr ,C 



Lit) 
L(t w ) 



(10) 



Here, following Fisher and Huse p!o| |, we assume 
the scaling function C(x) is related to a probability 
P s (L(i)/L(i w )) that a given spin belongs to the same 
domain at the two different epochs characterized by Lit) 
and L(i w ). It satisfies (7(1) = 1 and C(x) ~ x~ x at 
x ^> 1 with A being in the range d/2 < A < d. It should 
be emphasized that, in contrast to usual domain growth 
process jlgj , we put the dynamical order parameter 
rather than (7ea as the amplitude of (7(i, t w ) since we 
obtained C(t,t w )\ xr ^i ~ qu above. 

The ZFC susceptibility in the aging regime becomes, 



Xzfc(Mw) ~ x( L (T~),L(t)) 



Atw) L l ( L / L o) \ L J 



where c" is a numerical constant. The first term is due to 
the response of the droplets equilibrated within the tem- 
poral domain of size L(i) (see (Q)). The second term rep- 
resents contributions from the magnetizations (per spin) 
of order /5(0)/iTO 2 /T(L/L ) e due to droplets as large as 
L which are induced by the field in the past epochs, say 
around t with L(t) = L, and are depolarized by fur- 
ther domain growth. We also note that this term but 
with the integration from Lq to Lit w ) yields the TRM 
susceptibility written as JhJ xtrm(£ = t + t w ,t w ) ~ 
c nst piOW/T{Lit w )/L )- 9 {Lit)/Lit w ))- x with c nst = 
c" f dyy~ 1 ~ e+x > 0. These susceptibilities satisfy the 
sum rule xzfc(Mw) + Xtrm(Mw) = Xzfc(^O) which 
must hold for linear responses. 

In the special limit of t w = we obtain, 



XzFc(i,0) ~ xd - c" 



p(0)m 2 



T(L(t)/L ) e 



(12) 



with c'" — Ail) — Cnst. Note that xzfc(^O) converges 
to the dynamical susceptibility xd(> Xea) in the limit 
t — > oo. This implies xd is nothing but the field cooled 
(FC) susceptibility xfc @ i while xea is close to what is 
called the ZFC susceptibility xzfc- Thus we expect the 
well known experimental observation fi |) xfc > Xzfc 
is a dynamical but not a transient phenomenon . 

The overall feature of the two time quantities so far 
discussed is displayed in Fig. 111. Note that a paramet- 
ric plot of (7 vs xzfc becomes radically different from 
the conventional picture ^ in which the breaking of 
FDT and TTI are supposed to happen simultaneously at 
(<7ea,&bTxea)- 

Finally let us briefly discuss some numerical results 
concerning our scenario. We performed Monte Carlo 
(MC) simulations of isothermal aging of a 4 dimensional 
EA spin-glass model (size 24 4 )with ± J interactions (T c = 
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FIG. 2: ZFC susceptibilities and spin autocorrelation func- 
tions vs 1/L°(t). The symbols are (T/J)xzfc(t + t w ,t„) 
with t w = 0, 10, 10 2 , 10 3 , 10 4 , 10 5 MC step from the top to the 
bottom. The maximum time separation is r = 10 5 MC step. 
The curves with lines are corresponding 1 — C(r +t w , t w ) . The 
data of (T/J)xoq(i") (filled triangle) is shown at the bottom. 



2.0J, with k-Q = 1 here and hereafter) starting from ran- 
dom initial conditions. A set of data of 1 — C(r + t w ,t w ) 
and (T/J)xzFc(j~+tw,t w ) at T/J = 0.8 is plotted against 
L(t)-° in Fig. |. The field of strength h/J = 0.1 was 
used and the linearity of the responses was checked by 
the sum rule, xzFc(i, tw) +Xtrm(A t w ) = Xzfc(*,0). We 
used 9 = 0.82 obtained by a defect free-energy analysis 
p5[ . For the growth law L(t), we used the result of an in- 
dependent measurement of the growth of domain |2l], ^2) . 
As explained in jl4|, |T[ |2^] the analysis of the quasi- 
cquilibrium yields the equilibrium limit curve shown at 
the bottom of the figure (See It follows the scaling 

form @ pointing toward (T/J)xea — 0.18. The other 
extreme Xzfc( t i0) a ls° becomes linear in this plot as 
expected in ( |l2| ) pointing toward (T/J)xd — 0.34 well 
above xea- Apparently FDT ((T/J)x = 1 - C) is well 
satisfied at small r and broken at large r. Interestingly 
enough, the break points of the FDT move further away 
from (T/ J)xeq(r) suggesting the separation of the break- 
ing of FDT and TTI as we anticipated. Such a feature 
has not been reported in previous numerical studies [^6j . 

We also confirmed most of other scaling ansatz of the 
two time quantities within the 4D EA model which will 
be presented elsewhere [22[. It will be certainly inter- 
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esting to test our scenario by experiments such as noise 
measurements ]2^[ and further numerical simulations. 
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